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ABSTRACT 

We analyze the X-ray properties of a sample of local and high redshift galaxy clusters ex- 
tracted from a large cosmological hydrodynamical simulation. This simulation has been re- 
Q-f alized using the Tree+SPH code GADGET -2 for a ACDM model. It includes radiative cool- 

q ing, star formation and supernova feedback and allows to resolve radially the thermodynamic 

$_i ■ structure of clusters up to redshift 1 in a way that is not yet completely accessible to observa- 

£^ \ tions. We consider only objects with T cw > 2 keV to avoid the large scatter in the physical 

C3 ■ properties present at the scale of groups and compare their properties to recent observational 

constraints. In our analysis, we adopt an approach that mimics observations, associating with 
each measurement an error comparable with recent observations and providing best-fit results 
■ via robust techniques. Within the clusters, baryons are distributed among (i) a cold neutral 

phase, with a relative contribution that increases from less than 1 to 3 per cent at higher red- 
shift, (ii) stars which contribute with about 20 per cent and (iii) the X-ray emitting plasma that 
contributes by 80 (76) per cent at z = (1) to the total baryonic budget. A depletion of the 
cosmic baryon fraction of ^7 (at z = 0) and 5 (at z = 1) per cent is measured at the virial 
radius, i? v j r , in good agreement with adiabatic hydrodynamical simulations. We confirm that, 
also at redshift > 0.5, power-law relations hold between gas temperature, T, bolometric lumi- 
nosity, L, central entropy, S, gas mass, M gas , and total gravitating mass, Af to t and that these 
relations are steeper than predicted by simple gravitational collapse. A significant, negative 
evolution in the L~T and L — M to t relations and positive evolution in the S — T relation are 
detected at 0.5 < z < 1 in this set of simulated galaxy clusters. This is partially consistent 
with recent analyses of the observed properties of z > 0.5 X-ray galaxy clusters. By fixing 
the slope to the values predicted by simple gravitational collapse, we measure at high redshift 
normalizations lower by 10^4-0 per cent in the L — T, M to t — T, M gas — T, / gas — T and 
L — Aftot relations than the observed estimates. This suggests that either the amount of hot 
X-ray emitting plasma measured in the central regions of simulated systems is smaller than 
the observed one or a systematic higher value than actually measured of gas temperatures and 
total masses is recovered in the present simulated dataset. 

Key words: cosmology: miscellaneous - methods: numerical - galaxies: cluster: general - 
X-ray: galaxies. 



1 INTRODUCTION 

Only in the recent years, and after deep exposures obtained through 
the new veneration of X-rav relescnnes like Chandra and XMM- 
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Newton , it has been possible to start investigating systematically 
the physical properties of galaxy clusters, at redshift higher than 0.4 
(e.g. Holden et al. 2002, Vikhlinin et al. 2002, Ettori et al. 2004). 

Galaxy clusters are believed to form under the action of the 
dark matter gravity in the hierarchical scenario of cosmic structure 
formation. They assemble cosmic baryons from the field and heat 
them up through adiabatic compression and shocks that take place 
during the dark matter halo collapse and accretion. Simple self- 
similar relations between the physical properties in clusters are then 
predicted (e.g. Kaiser 1986, 1991, Evrard & Henry 1991) since the 
gravity does not have any preferred scale and the hydrostatic equi- 
librium between X-ray (mostly originated through bremsstrahlung 
processes) emitting gas and the cluster potential is a reasonable as- 
sumption. Therefore, how the baryons distribute within the cluster 
potential as function of cosmological time and whether or not scal- 
ing relations between integrated physical properties still hold and, 
possibly, evolve at high redshift, these are fundamental diagnostics 
to validate the self-similar scenario and to study the formation and 
evolution of X-ray galaxy clusters. A main role in assessing these 
issues is played by cosmological hydrodynamical simulations that 
can drive the analysis of observationally unresolved high redshift 
systems. 

To this purpose, we have analysed the predictions obtained at 
high redshift from the large hydrodynamical cosmological simula- 
tion presented and discussed in Borgani et al. (2004, hereafter Pa- 
per I), that has been shown to reproduce at z = several observed 
properties of X-ray clusters, such as the local temperature function, 
the luminosity-temperature relation for systems hotter than 2 keV, 
and the mass-temperature relation. 

As a natural follow-up of that work, we study here the evolu- 
tion of the X-ray properties of the simulated galaxy clusters from 
z = out to z = 1. This paper is organized as follows: in the 
next section, we describe our dataset of simulated clusters and their 
physical properties; in Section 3 we discuss how the baryons, and 
particularly the X-ray emitting plasma, are distributed in the local 
and z ^ 0.5 systems. In Section 4, we introduce the X-ray scaling 
relations and discuss their evolution at redshift 0.5 — 1 with a direct 
comparison with the observational results. Finally, we summarize 
and discuss our results in Section 5. 



2 PROPERTIES OF THE SIMULATED CLUSTERS 

The simulated clusters are extracted from the large-scale hydro- 
dynamical simulations described in Paper I. We refer to that pa- 
per for a detailed description of the simulation, while we provide 
here only a short summary. The simulated cosmological model is 
a standard flat ACDM universe, with fi m = 1 — Qa = 0.3, 
a 8 = 0.8, Q h h 2 = 0.019 and H = lOO/i km s" 1 Mpc~ x with 
h = 0.7. The simulation follows the evolution of 480 3 dark mat- 
ters and an initially equal number of gas particles within a box 
of 192/i~ 1 Mpc on a side, so that thdm = 6.6 x 10 9 M Q and 
m sss = 9.9 x 10 s Mq for the mass of the DM and gas particles, re- 
spectively. The Plummer-equivalent gravitational softening of the 
simulation was set to epi = 7.5 h" 1 kpc comoving from z = 2 
to z = 0, while it was taken to be fixed in physical units at higher 
redshift. 

The run has been realized using GADGET-2 1 , a massively par- 
allel tree N-body/SPH code (Springel, Yoshida & White 2001) 



1 |http : / /www .mpa-garching.mpg.de/gadgetl 
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Figure 1. Comparison between the true cluster masses and estimates mea- 
sured by using an isothermal /3— model. Mean deviations of about 15 per 
cent (at z = 0) or less (at z = 1) are obtained. 



with fully adaptive time-step integration, which includes the ex- 
plicit energy and entropy conserving SPH formulation by Springel 
& Hernquist (2002). The simulation includes radiative cooling 
of a plasma of primordial composition, the effect of a photoion- 
izing, time-dependent, uniform UV background (e.g., Haardt & 
Madau 1999). Star formation is treated using the self-regulated 
hybrid multiphase model for the interstellar medium introduced 
by Springel & Hernquist (2003). The code also includes a phe- 
nomenological description of galactic winds powered by type-II 
supernovae, such that star-forming gas particles contribute to the 
wind with a mass outflow rate two times larger than their star for- 
mation rate, with a wind velocity of about 360 km s _1 . 

Clusters are identified in the simulation box by first apply- 
ing a friends-of-friends halo finder to the distribution of DM par- 
ticles, with a linking length equal to 0.15 times their mean sep- 
aration. For each group of linked particles, we identify the parti- 
cle having the minimum value of the gravitational potential. This 
particle is then used as the center of the cluster to run a spheri- 
cal overdensity algorithm, which determines the radius around the 
target particle that encompasses a given overdensity. In the follow- 
ing analysis, we will consider a typical overdensity with respect to 
the critical density estimated at redshift z, p c ,z = 3H%/(8irG), of 
A z = 500 x A c , 2 /A e ,BdS; where A CiZ is the ratio between virial 
and critical overdensity at redshift z for our cosmological model, 
while A c ,Eds is the same quantity in an Einstein-de-Sitter cosmol- 
ogy. This overdensity is chosen for convenience, since the regions 
within these overdensities are the most well-studied ones from an 
observational point of view (e.g. Ettori et al. 2004). Hereafter, we 
indicate by "A = 500" the cluster regions enclosed within a sphere 
with radius Rsoo and total gravitating mass 



4 3 
M 50 o = -7rp c , z A z 7i5 00 . 



(I) 



We measure, on average, R500 to be 0.63 times i? v i r and constant 
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Figure 2. All (lines) and average (dots) differential 3-D gas temperature (upper panels) and density (lower panels) profile at redshift (panels on the left) and 
1 (panels on the right). The gas density is normalized to the virial value, p v i r , that is denned as the mean density within R v i r and is equal to p c , z X A z . 



in redshift (we define R v i r as the radius encompassing an average 
density equal to the virial density for our cosmology and at a given 
redshift; e.g., Eke et al. 1996, Bryan & Norman 1998). This con- 
firms that a fixed fraction of the virial radius is mapped by using the 
approximation to the solution for the collapse of a spherical top-hat 
perturbation. It is worth noticing that if one fixes the overdensity 
independently of the cosmological parameters and redshift, larger 
portion of the virial regions are mapped at higher redshift, since 
dark matter halos are more concentrated (i.e. have lower radial ra- 
dius). E.g., assuming A z = 500, one obtains Rsoo/Rvir = 0.49, 
0.56 and 0.58 at redshift 0, 0.5 and 1, respectively. 



As we show in Paper I and discuss further below in this sec- 
tion, the simulated gas temperature and density profiles present 
clear deviations from the mean behaviour of the observed ones. 
Therefore, we decide to adopt the mass as defined in equation Q 
and inferred from the simulated dark matter halos, instead of the 
corresponding value measured, e.g., through the hydrostatic equi- 
librium. We have measured, however, the reconstructed mass es- 
timates at 7?5oo by using an isothermal f3— model (Cavaliere & 
Fusco-Femiano 1976) with a single emission-weighted tempera- 
ture and inferring 7?soo from the corresponding mass profile (this 
procedure follows very closely what is generally done with obser- 
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Figure 3. Gas density profiles at high redshift as function of r/i?5oo- 
Dashed lines: best-fit profiles from the Chandra exposures of the 1 1 galaxy 
clusters with z > 0.8 (Ettori et al. 2004). Shaded region: (— 1, +1)<t range 
around the mean value of the 1 1 observed profiles. Dots: average (plotted 
here with the corresponding dispersion) of the simulated density profiles at 
redshift 1 (see the right lower panel in Fig.l2*l. 



Figure 4. Distribution of the different components of the baryonic budget 
within the cluster virial radius as a function of the total mass. The three 
concentrations of points correspond to (from the bottom to the top) / nou t 
(~2 per cent), / ncut + /star = /cold (about 20 per cent) and / cold + / gas 
(95 per cent) in unit of the cosmic baryon value Q^/Qm. Each component 
fi is equal to M;/M v i r . The plus symbols are for the z = objects, the 
squares for z = 0.5, the triangles for z = 0.7 and the crosses for z = 1. 



Table 1. Distribution of the baryons in the selected sample of simulated 
galaxy clusters. The mean values of the indicated ratios are quoted. 



n.obj 



97 0.007 
0.5 83 0.015 
0.7 81 0.021 

1 72 0.031 



0.194 
0.198 
0.201 
0.207 



(fW(n m ) 



0.799 0.930 

0.787 0.944 

0.779 0.949 

0.761 0.944 



vational data, e.g. Ettori et al. 2004). It is worth noticing that the 
mass profile is then a direct function of the best-fit modelization 
of the gas density profile that, as we discuss at the end of this sec- 
tion, shows evident mismatches with what is typically observed. 
We obtain that the /?— model estimates underestimate the true clus- 
ter mass by about 15 per cent on average at z — and are in good 
agreement at z = 1 even though with a larger dispersion around the 
mean value of the distribution of the relative deviations (see Fig.0. 
More relevant to the aim of this work is to note that these different 
definitions of the total gravitating mass, Af to t, do not affect sig- 
nificantly our results on the behaviour of the scaling relations. We 
refer to a forthcoming work in which the mass estimators applied 
to simulated and observed data will be discussed in detail (see also 
Rasia et al. 2004). 

According to the multi-phase model, each gas particle of suf- 
ficiently high overdensity is assumed to be composed of a hot ion- 
ized phase and of a cold neutral phase, whose relative amounts de- 
pend on the local conditions of density and temperature (Springel 



& Hernquist 2003). By its nature, the neutral cold component is 
assumed not to emit any X-rays. Following Croft et al. (2001), 
we also exclude those particles from the computation of X-ray 
emissivity having temperature below 3 x 10 4 K and gas density 
> 500/9 D ar, being p Da r the mean baryon density. In the following, 
we refer to cold gas as the gas which does not contribute to the X- 
ray emission and we distinguish it between a gas in a cold neutral 
phase, neut, and stars. 

The X-ray luminosity of each cluster in a given energy band 
is defined as 



L x = (fJ-nip) 



2 E 



m,h,iph,ih.{Ei, E 2 ; Ti] 



(2) 



where A(E\, E 2 ;T) is the cooling function in the energy band 
[EijEz]. In this equation the sum runs over all 7V ga s gas parti- 
cles falling within the cluster selected region, and fj, is the mean 
molecular weight (= 0.6 for a gas of primordial composition), m p 
is the proton mass, rrih.i and ph,% are the mass and the density asso- 
ciated with the hot phase of the i-th gas particle, respectively. The 
cooling function is computed from a Raymond-Smith code (Ray- 
mond & Smith 1977) by assuming zero metallicity. For the purpose 
of comparing with observational results on the scaling relations, in 
the following luminosities are always given in the bolometric band. 
We define the emission-weighted temperature, T e w, as 



J2f sa ° m Ki p Ki k(E x ,E 2 \ Tj) Tj 
£f gaB m h , i p h , i A(E 1 ,E 2 ;Ti) 



(3) 



We take Ei = 0.5(1 + z) keV and E 2 = 10(1 + z) keV so as to 
reproduce the observational procedure in the estimate of the tem- 
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perature from the measured photon spectrum, whose reconstruction 
at low energies, say below 0.5 keV, is made hard by instrumental 
limitations. Notice however that, using simulated Chandra obser- 
vations of galaxy clusters obtained with the software package X- 
MAS (X-ray MAp Simulator), Gardini et al. (2004) showed that 
the emission-weighted temperature inferred from hydrodynamical 
simulations can be significantly higher than the spectroscopic value 
extracted from observations (see also Mathiesen & Evrard 2001; 
Mazzotta et al. 2004). 

Finally, we define the entropy of the i-th gas particle as 
Si = Ti/n,i te 2 ^ 3 where n^e is the number density of free elec- 
trons associated to that gas particle. 

Out of the ensemble of simulated bound structures, we anal- 
yse the 97, 83, 81 and 72 systems at redshift 0, 0.5, 0.7 and 1, 
respectively, that have an emission weighted gas temperature, T ew , 
within A = 500 larger than 2 keV. This selection allows to con- 
sider only the objects where the gravitational collapse dominates 
the energetic budget and that present a local temperature function 
and luminosity-temperature relation in good agreement with the 
observed ones (see Paper I). 

In Figure|2| we plot the gas temperature and density profiles 
of the examined clusters at redshift and 1. At higher redshift, (i) 
the radial profiles tend to be more spread around the mean value, 
suggesting that many systems are still in formation, (ii) the tem- 
perature distribution appears slightly flatter within the virial radius, 
decreasing by about 50 per cent from 0.2 x i? v i r to 0.8 x R vlI 
(to be compared with a reduction by ~60 per cent at z = 0), 
and (iii) the central (r < 0.2 x R v i T ) gas density is less peaked 
by about 50 per cent than what observed at z = 0. However, the 
mean polytropic index 7, measured as the coefficient of the linear 
fit logTg a s(r) = k + (7 — 1) x logpg as (r), remains almost con- 
stant with redshift, varying from 1.18 (rms: 0.04) at z — to 1.16 
(rms: 0.06) at z = 1. We conclude that the overall shape of the 
gas temperature and density profiles does not change significantly 
with redshift, even though local radial variations are present as a 
consequence of the different dynamical state of the objects, with 
dynamically younger systems being located at higher redshift. 

In Figure [3] we show the comparison between simulated and 
observed gas density profiles at 0.8 < z < 1.3. The observed 
profiles are recovered from the best-fit results obtained by apply- 
ing a /3— model to deep Chandra exposures (Ettori et al. 2004) and 
are normalized by the corresponding p v i r (see caption in Fig. [2} 
as done for simulated profiles. The overall shape seems in good 
agreement with the mean estimated values, even though the scatter 
in the observed data (indicated by the gray region in Fig.[3j is defi- 
nitely larger than what is measured in the simulated ones at z = 1. 
This suggests that either the real data at high— z are subjected to a 
more complicated dynamical history than actually simulated or the 
latter ones are more regular due to the criteria adopted in select- 
ing the dark matter halos. In detail, the more significant deviations 
are within 0.1 x -R500 and above 0.7 x R500, where the steeper 
simulated profiles have mean values higher by 10-20 per cent and 
lower by 20 per cent (by 25 per cent at -R500) than the observed 
ones, respectively. The difference remains within few per cent in 
the range (0.1 — 0.7) x Rsoo- On the other hand, once we limit our 
analysis to the hottest systems by selecting the 10 objects with the 
highest T cw (< -R500), the discrepancy between the average pro- 
files becomes more noticeable: deviations larger than 10 per cent 
are present within 0.4 x -R500 and beyond 0.8 x Tisoo- However, 
once we compare the rough estimate of the expected mean gas mass 
by integrating these average profiles over the cluster volume up to 
-R500, we obtain a reasonable agreement between the observed and 



simulated values when only the hottest systems are analysed. A 
lower estimate by about 10 per cent is measured when the average 
gas density profile is obtained from all the simulated z — 1 objects. 

The few observational data available on the gas temperature 
profiles at high— z (e.g. Jeltema et al. 2001, Arnaud et al. 2002) 
indicate isothermality up to half the virial radius but with typical 
relative uncertainties of 20 per cent or more (90 per cent confidence 
level) on each temperature measurement of the 3-4 radial bins. Mis- 
matches between observed and simulated temperature profiles of 
the same order as highlighted in Paper I for local systems are prob- 
ably present here but can not be confirmed. We refer the reader to 
section 3.3 and 3.4 of that paper for a detailed discussion of these 
discrepancies. 



3 THE DISTRIBUTION OF BARYONS 

The distribution of the baryons, and of the hot X-ray emitting gas 
in particular, is a diagnostic of the energetic phenomena that take 
place during the formation and evolution of galaxy clusters. Since 
this distribution affects the estimates of the gas density, entropy and 
emissivity, the analysis of how baryons populate the cluster poten- 
tial well provides information on the properties of the integrated 
measurements. 

Overall, we find between 1 (at redshift 0) and 3 (at z = 1) 
per cent of the baryons to be present in the cold phase, about 20 
per cent in stars and between 77 (at 2 = 1) and 80 (at z — 0) 
per cent in the hot, X-ray emitting plasma. We also note a slight 
decrease from 7 to 5 per cent in the depletion of the baryons with 
respect to the assumed cosmic baryonic budget from z = to 1 (see 
Table 0. These values of depletion, that have to be considered as 
lower limits considering the measured over-production of stars, are 
however comparable to what predicted from adiabatic simulations 
such as the ones considered in the Santa Barbara Project (Frenk et 
al. 1999) that indicate a mean fb/(Slb/tom) of 0.92 (rms: 0.06). 

The overall amount of cold gas, Af ncut + M s tar is mainly 
composed of stars formed in high density regions and accounts for 
about 20 per cent of the total baryon budget (see average value in 
our samples in Table 0, at variance with a ~ 10 per cent esti- 
mate observed in nearby clusters (e.g. Lin et al. 2003). We point 
out that the major part of the star overproduction in simulated clus- 
ters takes place in the central cD galaxy, where the action of the 
galactic winds is not efficient enough to prevent overcooling. The 
lack of observational evidence for low-temperature gas in central 
cluster regions (e.g., Peterson et al. 2003) calls for the need of a 
central heating mechanism. We argue that this mechanism, which is 
not included in the present simulations, should reduce the resulting 
amount of stars. Moreover, there is a mild evidence for a higher effi- 
ciency to form stars in smaller systems (a trend also found in obser- 
vations, see Fig. 7 in Lin et al. 2003) and at higher redshift (Fig.|4}- 
As Springel & Hernquist (2003) suggest, this is an evidence for a 
less efficient formation of cooling flows in halos with virial temper- 
ature above 10 7 K. However, an analysis of our simulated clusters 
(Murante et al. 2004) shows that we have a diffuse star component, 
not associated with any galaxy, whose mass ranges from 10 per 
cent to more than 40 per cent of the total stellar mass of the cluster, 
being higher for more massive objects. Thus, the disagreement be- 
tween the amount of cold baryons present in simulations and those 
observed could be lessened. Furthermore, depending on the spa- 
tial distribution of the diffuse component, the estimates of the star 
formation efficiencies could be slightly biased both in observations 
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Figure 5. The gas fraction as a function of the radius in simulated clusters at z = (left) and z = 1 (right). The dots represent the best-fit obtained by using 
the functional form in equationRI 



and in simulations, in the direction of having a smaller efficiency 
for massive clusters. 

The X-ray emitting gas fraction, / gas (< r) = M gas (< 
r)/M to t(< r), increases when moving outward to (0.3 — 0.4) x 
.Rvir and flattens then toward a constant value that is about 0.8 x 
/bar ~ 0.75 x Ob/fim (Fig-EJ. This confirms that the hot ICM is 
more extended than the dark matter distribution, as already known 
from observational data (e.g. David et al. 1995, Ettori & Fabian 
1999), and that one needs to sample the cluster X-ray emission out 
to ~ 0.4 x 7?vi r « 0.7 x 7?5oo to fully recover the baryon fraction 
in the form of hot plasma. 

We find that a simple functional form can reproduce this de- 
pendence upon the radius quite well in the regions beyond the in- 
ner steep gradient, / gas (< r) oc (2r) v / (R v i T + r) v , with r\ about 
0.2-0.3, whereas a further suppression is required to reproduce the 
profile within ~ 0.3 x 7? v i r : 



/gas(< r) 
"/gas ( Rvir 



-.1 



2 r/R viv 
1 + r/R vil 



(4) 



with r s ~ 0.3-Rvh-- The best-fit results obtained by applying this 
functional form to the average radial profile of / gas are: / gas (< 
Rvii) = 0.75 x (f^b/fim), V between 0.24 and 0.32, and r s increas- 
ing from 0.28-Rvir to 0.34i? v i r when systems at higher redshift are 
considered. 

To explore the gas mass fraction as a function of radius in 
the hottest systems that are more accessible to the observational 
analysis, we select the 10 hottest systems at redshift and 1 (av- 
erage value of T cw (< R500) =4.8 and 3.4 keV, respectively). We 
measure then a mean (standard deviation) / gas / (^b/^m) of 0.638 
(0.030) and 0.567 (0.071) at ifesoo ~ 0.28/?^, 0.732 (0.020) and 
0.698 (0.026) at i? 5 oo, 0.760 (0.021) and 0.731 (0.015) within R vil . 
Even though a larger contribution of hot gas to the cosmic fraction 
is measured at higher redshift, these measurements are in reason- 
able agreement once the associated scatter is considered. 
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Figure 6. The gas fraction within A = 500 as a function of the gas temper- 
ature. The symbols of the points are as in caption of Fig.0 The horizontal 
lines indicate the mean values of the gas fraction at different redshift (dotted 
line: z = 0; dash-dotted line: z = 0.5; dashed line: z = 0.7; long dashes: 
z = 1). 



4 THE SCALING RELATIONS AND THEIR EVOLUTION 

Under the assumptions that the smooth and spherically symmet- 
ric intra-cluster medium (ICM) emits by bremsstrahlung and is 
in hydrostatic equilibrium with the underlying dark matter po- 
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Table 2. Best-fit results on the local (z 



0) scaling relations (dashed and dotted lines in the plots; see equation |5) of simulated galaxy clusters. The 

n44 „ rrr c — 1. t U~ tr.tr.} m-acc A/f. . in 1 fll4 U„ . fh^ nnc mace A if 4n 1 0.13 J 



temperature, T, is in unit of 6 keV; the bolometric luminosity, L, in 10 44 erg s 1 ; the total mass, A/tot, in 10 14 M Q ;the gas mass, Aig as , in 10^ Mq. All 
these quantities are estimated within Rsoo- The entropy, S, is in unit of keV cm -2 and is measured at 0.1 X i?200- When the slope A is fixed, we estimate 
the eiTor-weighted mean of (log Y — A log X ) and evaluate the error after resampling Y and X by 1 000 times according to their uncertainties. The best-fit 
on the evolution parameter B (see equation[6| is obtained by considering (i) all the objects at redshift 0.5, (ii) those at ^ 0.7 and (iii) only the clusters at 
2 = 1. The quoted errors are at la level (68.3 per cent level of confidence and A\ 2 = 1 for one interesting parameter). 



relation (Y - X) 



B(z^ 0.5) B (z > 0.7) B (z = 1) 



F z M t , 



1.07(±0.06) 
0.87(±0.01) 



2.08(±0.15) 
1.50 (fixed) 



-0.12(±0.05) 
-0.20(±0.04) 



-0.14(±0.05) 
-0.22(±0.04) 



-0.20(±0.07) 
-0.30(±0.05) 



F Z Mr, 



1.14(±0.08) 
0.83(±0.01) 



2.40(±0.20) 
1.50 (fixed) 



-0.12(±0.05) 
-0.26(±0.03) 



-0.14(±0.06) 
-0.30(±0.03) 



-0.22(±0.08) 
-0.38(±0.04) 



1.23(±0.14) 
0.76(±0.01) 



3.33(±0.35) 
2.00 (fixed) 



-0.76(±0.08) 
-0.96(±0.04) 



-0.76(±0.09) 
-0.98(±0.04) 



-0.80(±0.11) 
-1.04(±0.05) 



F Z M U 



-0.47(±0.04) 
-0.39(±0.02) 



1.57(±0.11) 
1.33 (fixed) 



-0.56(±0.05) 
-0.68(±0.04) 



-0.56(±0.05) 
-0.66(±0.04) 



-0.50(±0.07) 
-0.64(±0.05) 



-0.35(±0.03) 
-0.34(±0.01) 



1.38(±0.07) 
1.33 (fixed) 



-0.58(±0.03) 
-0.60(±0.02) 



-0.56(±0.03) 
-0.58(±0.03) 



-0.50(±0.04) 
-0.54(±0.04) 



-1.01(±0.01) 
-1.04(±0.01) 



0.09(±0.02) 
0.00 (fixed) 



-0.04(±0.03) 
-0.06(±0.03) 



-0.06(±0.04) 
-0.06(±0.04) 



-0.08(±0.05) 
-0.08(±0.04) 



? 4 /3< 



2.81(±0.06) 
2.66(±0.01) 



1.42(±0.17) 
1.00 (fixed) 



0.26(±0.04) 
0.22(±0.02) 



0.24(±0.04) 
0.18(±0.02) 



0.18(±0.05) 
0.12(±0.03) 



tential, one can relate the bolometric luminosity, L, gas temper- 
ature, T, gas entropy, S, gas mass, M gas , and total mass, M to t 
in a simple and straightforward way. For instance, the equation 
of hydrostatic equilibrium, d(p gas T)/dr ~ p gas GMt t/r 2 , al- 
lows to write M tot oc TR that, once combined with the defi- 
nition of the total mass within a given overdensity A z , Mtot oc 
E 2 Z A Z R 3 , implies E z &l /2 M tot oc T 3/2 , where E z = H z /H = 



[n m (l + z) A + 1 — f2 m ] (for a flat cosmology with matter 
density fi m and Hubble constant at the present time Ho). Similarly, 
the definition of the bremsstrahlung emissivity e oc A(T)ng aa oc 



T l/2 r 



(the latter being valid for systems sufficiently hot, 



e.g. > 2 keV) allows to relate the bolometric luminosity, L, 



and the gas temperature, T: L 

T 1/2 f*M? oi R- 



eR J 



fga. s T , where we have made use of the for- 
mula on the Mtot — T relation shown above. 

By combining these basic equations, we can obtain the scaling 
relations among the X-ray properties that we are going to investi- 
gate in the present section: 

• F z Mtot oc T z ' 2 
. F z M gas oc T 3 / 2 

. / gas oc T° oc Mt°ot 

• F~ z 

• f: 



L oc T 
S 1 L 

F z /3 S oc T, 



M tot ) 4/3 oc 



(F z M g . 



S4/3 



where we have combined all the cosmological dependence in the 
factor F z — E z x (A z /A z= o) 1 ^ 2 . All the quantities are measured 
within regions with overdensity A z = 500 (see Sect. 2), apart from 
the entropy S that is related to the thermodynamical entropy K ac- 



cording to K cc log S, and is estimated as T gas (i?o.i)/n 2 ^ 3 (7?o.i) 
withi?o.i = 0.1 x i?200- 

The relative error (at la level of confidence) on the physi- 
cal quantities is assumed to be 10 per cent for the luminosity, 15 
per cent for the temperature and gas mass and 25 per cent for the 
gravitational mass. These uncertainties are the average values cited 
presently for local and high— z clusters (e.g. Finoguenov, Reiprich 
& Bohringer 2001, Ettori et al. 2004). 

The behavior of these scaling laws are examined first in their 
normalization and slope by fitting the logarithmic relation 



log Y = a + A log X 



(5) 



between two sets of measured quantities {Xj} and {Y,}. We use 
the bisector modification (i.e. the best-fit results bisect those ob- 
tained from minimization in vertical and horizontal directions) of 
the linear regression algorithm in Akritas & Bershady (1996 and 
references therein, hereafter BCES for Bivariate Correlated Errors 
and intrinsic Scatter) that takes into account both any intrinsic 
scatter and errors on the two variables considered as symmetric. 
We adopt the best-fit results from the BCES(F|X) estimator only 
when the expected slope is zero, being the inverse one equals to in- 
finity, like in the cases of the / gas — T relation. The uncertainties on 
the best-fit results are obtained from 10000 bootstrap resampling. 

With a relevant difference with respect to the analysis pre- 
sented in Paper I, we have selected for the present analysis only 
galaxy clusters with T cw (J?5oo) > 2 keV. Furthermore, because 
of the inclusion of a relative error on the measured quantities, the 
best-fit results do not necessarily coincide with those provided in 
Paper I. We refer to the latter work for a detailed discussion of the 
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Table 3. Best-fit results on the simulated and observed high— z (z ^ 0.5) samples. The relations and the parameters are defined in Table |2| The evolution 
parameter B is estimated by using the local relations indicated in Tablel2lfor simulated objects, and the best-fit values in Finoguenov et al. (2001, Mtot — T), 



Mohr et al. (1999, M e . 



T and / gl 



T), Ettori et al. (2002, L - T), Reiprich & Bohringer (2002, L - M tot ) and Ponman et al. (2003, S - T) for 



the observed 22 clusters with z 0.5 (Ettori et al. 2004). The deviations in a between the parameter B as estimated in the simulations and observations are 
indicated in the column "erg". 



relation (Y - X) 



B 



from simulations 



B 



from observations 



0"B 



F z M t , 



F^L - T 
F Z ~ 1 L- F z M tot 



F~ X L - F z M gas -0.48(±0.01) 
-0.48(±0.01) 



F? /3 S - T 



1.29(±0.12) 2.72(±0.31) -0.12(±0.05) 0.85(±0.04) 1.89(±0.45) 0.24(±0.25) -1.44 

0.82(±0.01) 1.50 (fixed) -0.20(±0.04) 0.88(±0.02) 1.50 (fixed) 

1.49(±0.18) 3.35(±0.46) -0.12(±0.05) 0.83(±0.04) 2.63(±0.36) -0.52(±0.28) 1.42 

0.83(±0.01) 1.50 (fixed) -0.26(±0.03) 0.98(±0.02) 1.50 (fixed) 

1.45(±0.22) 4.33(±0.54) -0.76(±0.08) 0.58(±0.07) 4.12(±0.75) -2.20(±0.46) 3.12 

0.55(±0.01) 2.00 (fixed) -0.96(±0.04) 0.75(±0.03) 2.00 (fixed) 

-0.61(±0.02) 1.57(±0.07) -0.56(±0.05) -1.33(±0.60) 2.24(±0.71) -2.28(±0.31) 5.44 

-0.55(±0.01) 1.33 (fixed) -0.68(±0.04) -0.38(±0.04) 1.33 (fixed) 

1.31(±0.04) -0.58(±0.03) -0.70(±0.06) 1.55(±0.07) 

1.33 (fixed) -0.60(±0.02) -0.49(±0.02) 1.33 (fixed) 

-0.96(±0.03) 0.24(±0.09) -0.04(±0.03) -1.03(±0.04) 0.60(±0.49) -0.64(±0.22) 2.75 

-1.05(±0.01) 0.00 (fixed) -0.06(±0.03) -0.89(±0.02) 0.00 (fixed) 

3.19(±0.17) 2.22(±0.42) 0.26(±0.04) 2.75(±0.02) 0.70(±0.23) 0.68(±0.06) -5.87 

2.71(±0.01) 1.00 (fixed) 0.22(±0.02) 2.65(±0.01) 1.00 (fixed) 



(mis)matches between observed and simulated scaling relations at 
z = 0. We concentrate here on the properties of these relations at 
high redshift (z > 0.5). 

In the following sections, we discuss the slope, normalization 
and evolution with redshift of the power-law fit among these quan- 
tities and compare these results to the observational constraints ob- 
tained by Ettori et al. (2004) for a sub-sample of 22 X-ray galaxy 
clusters at z ^ 0.5 observed on kpc scales through deep Chandra 
observations. It is worth mentioning that, in the present analysis, we 
adopt systematically the definitions and methods discussed there to 
study the behaviour of the scaling laws of X-ray galaxy clusters. 
The results on the best-fit normalization and slope for the local scal- 
ing laws here investigated are given in Tableland shown as dotted 
(when the slope is fixed to the value predicted by the above scaling 
relations) and dashed (when the slope is a free parameter) lines in 
Figuresl7l- I12l In Table|3| a synoptic summary of the best-fit re- 
sults obtained from samples of simulated and observed clusters at 
z ^ 0.5 is given, with an assessment on the deviation between the 
evolution parameter B (see Section 4.2) as measured in the simu- 
lations and from observational data. 



4.1 On the slope and normalization of the scaling relations 

We fit equation[5]to the quantities of the simulated clusters at z = 
and quote the results in Table|2| Overall, we observe a slope steeper 
than predicted from self-similar model in the examined scaling 
laws, with L-T (Fig.0, Mtot - T (Fig.OQ and M gas - T (Fig.g) 
relations being the ones with deviations larger than 3<r (best-fit of 
A = 3.3 ± 0.3, 2.1 ± 0.2 and 2.4 ± 0.2, respectively), whereas 
S — T (Fig. ll2t deviates by about 2.5<r from the expected values 



of A = 1, L - Mtot by 2.2a and L 
predicted slope of 4/3. 

About the normalization of these scaling relations at high red- 
shift (z ^ 0.5) and how they compare to the observed estimates, we 
use the results by Ettori et al. (2004) once the slope is fixed to the 
expected value from simple gravitational collapse assumption. The 
simulated clusters present normalizations that are definitely lower 
than what is actually observed in the same redshift range (see Ta- 
ble|3j: ot is lower by 7 per cent in the Mtot — T relation, 12 per cent 
in the M gas — T relation, 27 per cent in the L — T relation, 45 per 
cent in the L — M to t relation, 18 per cent in the / gas — T relation and 
consistent within few per cent in the S — T relation. One possible 
explanation for these lower normalizations is that the amount of the 
hot X-ray emitting plasma measured within -R500 of simulated sys- 
tems at z ^ 0.5 is smaller than the observed one. This deficiency in 
the density of cosmic baryons detectable as hot gas in hydrodynam- 
ical simulations of galaxy clusters might be complementary to the 
well-known "over-cooling" problem, where radiative cooling con- 
verts a too large fraction of cosmic baryons into collisionless stars 
(e.g. Suginohara & Ostriker 1998, Balogh et al. 2001, Borgani et 
al. 2002) if the local cooling time is not conveniently increased by 
some sort of feedback process. On the other hand, the good agree- 
ment in the L — Af gas relation suggests that overestimates of the 
gas temperature and total mass with respect to what actually mea- 
sured can affect our overall results. We have indeed verified that by 
reducing, for example, by ~15 per cent the simulated estimates of 
both T and Mtot provides a general good match between the nor- 
malizations of the scaling relations in the observed and simulated 
datasets. We refer to a future work where we will compare differ- 
ent definitions of gas temperatures, as inferred from simulations, 
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Figure 8. The same as in Fig0 but for the Mtot — T relation. 
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Figure 7. The local L — T relation estimated at A = 500 as best-fit from 
the simulated clusters at z = (plus symbols). Dotted line: slope fixed 
to the predicted value. Dashed line: slope free. The simulated objects at 
z = 0.5 (squares), z = 0.7 (triangles) and z = 1 (crosses) are compared 
to the observed dataset from Ettori et al. (2004) with 0.5 < z < 0.8 (small 
dots) and z > 0.8 (large dots). (Lower panel) Ratio between the measured 
luminosities at z = (plus symbols) and z = 1 (crosses) and the best-fit 
local power-law. 
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Figure 9. The same as in Fig0 but for the Af gas — T relation. 
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Figure 10. The same as in FigQ but for the L — Mtot relation. 



to observational data, and the implications on scaling relations and 
measurements of the cluster mass. 



4.2 On the evolution of the scaling laws 

We measure a median gas temperature of 2.5, 2.4, 2.3 and 2.3 keV 
in our samples at redshift of 0, 0.5, 0.7 and 1, respectively. On 
the other hand, we note a drastic decrease in the median values of 
F Z _1 L/(10 44 erg s _1 ) from 1.01 at z = to 0.61, 0.55 and 0.42 at 
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Figure 11. The same as in Fig0 but for the L — Af gas relation. 
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Figure 12. The same as in FigQ but for the S — T relation evaluated at 
Ro.i = 0.1 x R 2 oo- 



z = 0.5, 0.7 and 1, mostly due to the cosmological correction in- 
troduced by the evolution of the Hubble constant and parametrized 
through E z . Significantly lower values than those measured locally 
are also evaluated for the quantity F z M to t / (1O 14 M0) that changes 
from 1.94 (z = 0) to 1.86 (z = 0.5), 1.63 (z = 0.7) and 1.53 
(z — 1). These simple calculations give already indications for a 
declining luminosity and mass for given temperature and with in- 
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Figure 13. Distribution of the values of the evolution parameter B as 
measured in sub-samples of 22 simulated objects to mimic the dataset at 
z 0.5 in Ettori et al. (2004). Dark shaded region: 1 — a range and 
best-fit (dashed line) of B as evaluated for the overall sample of simulated 
clusters. Light shaded region: 1 — a range and best-fit (dotted line) of B as 
evaluated in Ettori et al. (2004) for the same relation. 



creasing redshift or, in other words, for a negative evolution of the 
corresponding scaling relations. 

To constrain properly the evolution in the normalization of the 
considered scaling laws, we proceed as discussed in Ettori et al. 
(2004). First, we fix (a, A) to the best-fit results obtained from 
the clusters at z = 0, (a, A), and, then, evaluate the confidence 
interval through a least-square minimization on the parameter B in 
the relation 



logF = a + AlogX + Blog(l + z). 



(6) 



In particular, for a given grid of values of {-B;}, we search for the 
minimum of the merit function 



E 



[log Yj-a-A log Xj - Bi log(l + Zj)] ' 



(7) 



where the errors on the best-fit local values, e— and e-j, are 
considered and propagated with the uncertainties, ei og x = 
ex I (X In 10) and ei og y — ey/{Y In 10), on the measured quan- 
tities. 

We measure a small, but highly significant, negative (i.e. B < 
0) evolution in the L - T (B » -0.8) and L - Aftot (B « -0.6) 
relations. Marginally (< 3a) negative values of B are measured 
in the M tot ~ T, Af gas — T and / gas — T relations. A remark- 
able positive evolution is estimated in the S — T relation. All these 
trends becomes more significant when objects at higher redshift are 
considered. In particular, the relative differences between local and 
z = 1 objects have been plotted in the lower panels of FiguresQ- 
1121 where the corresponding physical quantities are referred to the 
best-fit local power-laws. A clear segregation between local and 
high— z clusters is present in, e.g., the L — T (Fig.|7J and L — M to t 
(Fig.[l0j plots. 
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All these values in the evolution of the studied scaling rela- 
tions are qualitatively consistent with the recent observational con- 
straints (e.g. Ettori et al. 2004 and see Table[3}. On the other side, 
the amount of measured evolution is less than actually observed in 
the simulated L — T, L — M to t and S — T relations. 

To check against any bias effect due to the number of objects 
selected in the sample, we have randomly extracted 7 objects at 
z = 0.5, 10 at z — 0.7 and 5 at z = 1 to mimic the distribu- 
tion of the 22 observed clusters at z ^ 0.5 in Ettori et al. (2004) 
in which 7 clusters have redshift between 0.5 and 0.65, 10 are be- 
tween z = 0.65 and z — 0.85 and 5 have z between 0.85 and 1.3. 
We have repeated the sampling for 1,000 times and plot the dis- 
tribution of the values of the evolution parameter B in fig. 1 131 We 
conclude that no-significant deviations from the best-fit B value 
as measured for the overall sample are obtained when sub-samples 
are considered. Moreover, the observed evolution in the L — T re- 
lation cannot be reproduced at very high significance (> 99.9 per 
cent), whereas the one observed in the Mtot — T relation is only 
marginally consistent with the distribution of the B values (about 
23 per cent of the simulated sub-samples have an evolution param- 
eter that falls within the 1 — a error range of the observed B value). 



4.3 On the results obtained by adopting different selection 
criteria 

We have repeated the analysis discussed above by extracting the 
simulated objects according to different criteria to answer two main 
issues: 

(i) . What happens if we select only very massive systems (e.g. 
Tew > 3 keV) in accordance to the observational datasets avail- 
able ? By construction, no very massive virialized systems are build 
in this cosmological simulations implying that a selection at higher 
threshold in temperature (mass) reduces significantly the number 
of objects available. When we select the simulated clusters with 
Tew > 3 keV, we collect just 24 objects at z = and 8 at z = 1 
(14 and 12 at z — 0.5 and 0.7, respectively). Therefore, by de- 
creasing the number of objects, we increase the statistical uncer- 
tainties of the best-fit parameters and no more significant evolution 
are detected. For example, in the L — T relation we constrain the 
local slope to be 3.4 ± 1.6, with a marginally negative evolution at 
z ~^ 0.5 (B ~ —0.7 ± 0.4). An evolution consistent with zero is 
measured in the M to t - T relation (B = -0.1 ± 0.3 at z ^ 0.5), 
which presents a slope at z — consistent with the predicted value 
of 1.5 (2.1 ± 1.1). 

(ii) . How do the relations change when the gas mass-weighted 
temperature, T mw , is used instead of the emission-weighted one 
(see, e.g., Allen et al. 2001, Thomas et al. 2002) ? If we select 
simulated clusters with T mw > 2 keV, we reduce significantly the 
number of objects selected (42 vs. 97 simulated clusters at z — 
and 12 vs. 72 at z = 1) and confirm by an indirect way the well- 
recognized problem of simulated cores being hotter than actually 
observed (see discussion on the temperature profile in section 3.4 
in Paper I). All the scaling relations appear to be flatter and more 
consistent with self-similar predictions than measured when T" e w is 
in use: e.g., slopes of 3.0 ± 0.6 and 1.6 ± 0.2 are measured in the 
local L — T and Mtot — T relations. Also the evolution parameter 
B departs less significantly from zero, with values that are typically 
50 per cent lower than the measured ones in the original selection. 



5 SUMMARY AND DISCUSSION 

Using a large cosmological hydrodynamical simulation, we have 
selected galaxy clusters at redshift z = 0, 0.5, 0.7 1, all having 
T cw within R500 larger than 2 keV, to explore the distribution of 
baryons and the X-ray scaling relations as function of redshift. We 
adopt an approach in our analysis that mimics observations, asso- 
ciating with each measurement an error comparable to recent ob- 
servations and providing best-fit results obtained with robust tech- 
niques. 

Our main findings can be summarized as follows: 

(i) . At z ^ 0.5, the simulated gas temperature and density pro- 
files tend to show more scatter around the mean value, suggesting 
that many systems are still in formation and dynamically younger 
than the ones simulated locally. The temperature and density gra- 
dients are flatter and less centrally peaked, respectively, than those 
at z = 0. However, the overall shape of these radial profiles repre- 
sented by the polytropic index does not change significantly. 

(ii) . Gas density profiles for simulated clusters are slightly 
steeper that actually observed at z ~ 1, with deviations of about 10- 
20 per cent (on average) below 0.1 x R500 and above 0.7 x i?5oo- 
This implies a mean integrated value of the simulated gas mass 
lower by about 10 per cent. A significantly larger scatter in the 
observed data suggests that either the complex dynamical history 
of the formation of these structures is not recovered fully in this 
hydrodynamical simulation or the simulated dark matter halos are 
more regular due to the criteria adopted in their selection. 

(iii) . The baryons within the virial radius are distributed among 
a cold phase, with a relative contribution that increases from less 
than 1 to 3 per cent at higher redshift, a collisionless phase in stars 
of about 20 per cent and the X-ray emitting plasma that contributes 
by 80 (76) per cent at z = 0(1) to the total baryonic budget. A 
depletion of cosmic baryons of the order of ~7 per cent is measured 
locally at the virial radius, R v i r . The average value of the depletion 
decreases slightly to 5 per cent at z — 1. These values are consistent 
with the results obtained in adiabatic hydrodynamical simulations, 
even though they have to be considered as lower limits considering 
the measured over-production of stars. 

(iv) . The X-ray emitting gas mass fraction, / gas (< r), increases 
with rasa power-law with a characteristic scale that also increases 
from 0.28 to 0.34 xi? v i r when going from z = to z = 1, with a 
typical slope of about 0.25. The measured gas fraction in the ten 
hottest systems at -R2500, ^500 and i? v ir is 0.62, 0.73 and 0.76 
x (n b /fi m ) at z = and 0.57, 0.70 and 0.73 x (n b /n m ) at * = 1, 
respectively. The differences as function of redshift are not signifi- 
cant being within the measured scatter. 

(v) . We confirm that, also at high— z, simulated clusters have X- 
ray scaling relations between temperature, T, luminosity, L, central 
entropy, S, gas mass, M gas , and total gravitating mass, Mtot, which 
are steeper than predicted by simple gravitational heating. 

(vi) . When we fix the slope of these scaling laws to the values 
expected from the self-similar scenario, normalizations a lower by 
10-40 per cent are measured in the L — T, Mtot — T, M gas — T and 
L — Mtot relations at z ^ 0.5 when compared to the observational 
results for high redshift clusters (e.g. Ettori et al. 2004). The nor- 
malization of the central entropy - temperature relation is instead 
higher by few per cent, confirming the lack of gas mass in the cen- 
tral regions of these high— z simulated systems. On the other hand, 
the good agreement in the L — M gas relation suggests that overes- 
timates of the gas temperature and total mass with respect to what 
actually measured can affect our overall results. We have indeed 
verified that by reducing by 15 per cent the simulated estimates of 
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both T and Mtot provides a general good match between the nor- 
malizations of the scaling relations in the observed and simulated 
clusters. In this perspective, it should be noticed that Mazzotta et 
al. (2004) have recently suggested a spectroscopic-like definition of 
the ICM temperature from simulations that better matches the tem- 
perature obtained from the spectral fit of observational data. This 
alternative definition tends to provide temperature estimates lower 
than the emission-weighted values when the thermal structure of 
the cluster in exam deviates more from an isothermal modelization. 

(vii). We measure a declining luminosity and mass for given 
temperature and with increasing redshift. This negative evolution 
is parametrized through a power-law dependence upon the redshift, 
(1 + z) B , with B w -0.8 for the L - T relation, B « -0.2 for 
the Mtot - T relation, B w -0.2 for the M gas - T relation and 
B ~ —0.6 for the L — M to t and L — M gaa relations. A slightly 
positive evolution of B ~ +0.2 is measured in the S — T relation. 
As shown in Table [3] these results are consistent (deviation less 
than 3a) with the observed evolution evaluated in a similar redshift 
range (e.g. Ettori et al. 2004) for the M to t - T, M gas - T and 
/gas — T relations. Significantly less evident evolution are instead 
measured in the simulated S — T, L— Mtot and L — T relations. We 
have verified that the amount of evolution estimated here is robust 
against bias effects due to the number of objects examined. 

Our results demonstrate that the present dataset of simulated 
galaxy clusters tends to evolve at higher redshift, with a flatter gas 
temperature profile and a less centrally peaked density profile and a 
statistically significant deficiency of X-ray emitting plasma (items 
i, iii, iv and vii above). 

The comparison with the observational constraints at z ^ 0.5, 
performed by adopting the same definitions of the examined quan- 
tities and analyzing comparable samples both in size and physi- 
cal properties, further indicates that the measured lower normal- 
izations (item vi above) suggest either a significant lack of X-ray 
luminous baryons or an overestimates of the simulated values of 
the gas temperatures and total gravitating mass. An evolution to- 
ward less luminous and massive systems for given gas temperature 
at higher redshift (item vii) is measured in the simulated dataset in 
good agreement with recent observational results. 

Overall, our results show how the study of the evolution in 
the scaling relations can give indications on the properties of the 
feedback that can be investigated through extensive hydrodynami- 
cal simulations. In particular, these results support the request for 
either more efficient feedback or extra physical processes, which 
are able to suppress (i) the over production of stars in the central 
cluster regions and (ii) the lack of cosmic baryons in the hot X- 
ray emitting phase. Moreover, this additional source of heating is 
required not only to provide extra energy but also to distribute it 
radially in such a way (iii) to improve the agreement between the 
simulated and observed gas temperature both locally and up to red- 
shift of 1, (iv) to reduce the steepening of the gas density profile at 
high— z, (v) to 'soften' the gas density profiles of low temperature 
objects by increasing the central entropy, as already discussed in 
Paper I. 

The work presented here demonstrates that we are approach- 
ing an age in which code efficiency and supercomputing capa- 
bilities make the simulations able to describe the formation and 
evolution of cosmic structures over a fairly large dynamical range 
in a fruitful symbiosis with the present observational limits. The 
real challenge for improving the exchange between numerical and 
observational cosmology is, on the one side, to build algorithms 
that more faithfully incorporate all those astrophysical processes 



that are suspected to give relevant contribution to the observational 
properties of highly non-linear and over-dense structures, and on 
the other side, to resolve spatially with a sufficient number of pho- 
tons even at z » 1 all the physical quantities of interest, starting 
with the temperature of the X-ray emitting plasma. 
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